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When a phase-separated binary (A + B) mixture is exposed to a wall, that preferentially attracts 
one of the components, interfaces between A-rich and B-rich domains in general meet the wall making 
a contact angle 0. Young's equation describes this angle in terms of a balance between the A — B 
interfacial tension jab and the surface tensions j w a, j w b between, respectively, the A- and B-rich 
phases and the wall, jab cos = j w a — JwB- By Monte Carlo simulations of bridges, formed by one 
of the components in a binary Lennard-Jones liquid, connecting the two walls of a nanoscopic slit 
pore, is estimated from the inclination of the interfaces, as a function of the wall- fluid interaction 
strength. The information on the surface tensions j w a, Jwb are obtained independently from a new 
thermodynamic integration method, while jab is found from the finite-size scaling analysis of the 
concentration distribution function. We show that Young's equation describes the contact angles 
of the actual nanoscale interfaces for this model rather accurately and location of the (first order) 
wetting transition is estimated. 

PACS numbers: 68.08Bc, 05.70.Np, 64.75.Jk 



Introduction and Motivation: Pure fluids and 
fluid mixtures confined to nanopores are of substantial 
recent research interest due to various applications, e.g., 
extraction of oil or gas from porous natural rocks, use as 
"molecular sieves" to separate fluids, and various appli- 
cations in nanofluidics such as "lab on a chip" -devices 
and preparation of materials in nanotechnology, etc. 0- 
At the same time, the interplay between finite-size 
effects, surface effects due to the pore walls, and in- 
terfacial phenomena pose challenging problems to the 
theoretical understanding on the basis of statistical me- 
chanics [7l4lo|. A fundamental concept in this context 
is the contact angle with which an interface between 
coexisting (bulk) phases in a phase-separated system 
meets a wall under conditions of incomplete wetting 0- 
13|- In the limit of macroscopically large droplets, for 
which the droplet radius is many orders of magnitude 
larger than any interfacial width, one believes that can 
be expressed in terms of a balance between appropriate 
interfacial tensions via Young's equation [14]. For the 
case of a binary (A + B) fluid mixture undergoing phase 
separation in the bulk into A-rich and 5-rich coexisting 
phases, this equation reads 



jab cos j wA - j u 



(1) 



where jab is the interfacial excess free energy due to an 
infinitely extended flat interface between the coexisting 
bulk phases. In Eq. (pQ), the excess free energies of the 
A-rich and B-rich phases due to the wall are denoted 
by J W A and j wB , respectively. 



While numerous techniques exist for precise exper- 
imental measurements of both the contact angle of 
macroscopic droplets and the interfacial tension jab, 
it is challenging to measure j w a, JwB independently 
with sufficient precision. Thus a test of Eq. (pQ) even 
for macroscopic droplets is a very difficult task and fur- 
thermore, often hampered by substrate inhomogeneities 
causing contact angle hysteresis [l], [TlJ • For pores in the 
nanoscale size range, the problem is much harder, since 
the atomistically diffuse character of the interfaces can 
no longer be neglected, and of course, in many circum- 
stances interfaces are curved where the understanding 
of how the interfacial tension jab(R) varies with the 
radius of curvature is still incomplete jl5l4l8| . An addi- 
tional complication already comes into play for sessile 
mesoscale droplets attached to walls: then Eq. (pQ) is 
modified by a line tension (r) contribution from the cir- 
cumference length 2nr of the sphere-cap shaped droplet 
right at the wall [H, [lo|, 



jab cos = j wA - JwB + r/r. 



(2) 



Estimation of the magnitude of the line tension r is 
another problem very controversially discussed in the 
literature (2ll-[25|. Indeed the proper identification of 
the line tension effects is subtle [26j. Despite a large 
experimental and theoretical activity [27]-[37| for study- 
ing capillary condensation, liquid bridges in slit pores, 
etc., typically the information obtained does not yield 
both the contact angle and all the excess free energies 
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in Eqs. (pQ) and (|2j) with the notable exception of recent 
simulation studies of the Ising lattice-gas model 3fl [37j . 



The purpose of the present work hence is to provide 
a first complete test of the validity of Eqs. (pQ) and ((2]) 
in the nano-scale via Monte Carlo simulations of an off- 
lattice model of a fluid binary mixture. We shall study 
the structure of liquid bridges of A-rich phases in a B- 
rich background for various choices of the interactions 
between the particles and the walls of a nanoscopic slit 
pore and extract effective contact angles 6 e R to describe 
the inclined interfaces that occur. In addition, we com- 
pute jab for bulk interfaces and develop a new method 
to compute 7 w a, JwB for pure semi-infinite phases ex- 
posed to corresponding walls. Extending the approach 
of Winter et al. we are also able to derive the 



line tension r as function of the contact angle 6. 

Model and Simulation Technique: Let us con- 
sider a binary (A + B) fluid consisting of N point parti- 
cles, labelled by index i at positions r$, in a box of linear 
dimensions L x L x D, with periodic boundary condi- 
tions in x- and y- directions and having impenetrable 
walls of area L x L each at z = and z = D. The par- 
ticles interact via pairwise potentials u(rij = \fl —fj\), 
constructed from the full Lennard- Jones (LJ) potential 
<t>LAriA = 4e a p[((Tap/rij) 12 - (^/r^-) 6 ], a, P G A, B, 
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u(rij < r c ) = 0Lj(rij)-(i)Lj(rc)-(rij-r c )- 



(3) 

while u{rij > r c ) = 0. This ensures that both the 
potential and the force are continuous at the cut-off 
distance r^- = r c . The potential parameters are chosen 
such that the mixture is symmetric: a a A — &bb — 
(Jab = cr, £aa = £bb = 2sab — £ and r c was set at 
2.5cr. Working at a reduced density p* = 1, where p* = 
pa 3 = Na 3 /V, neither crystallization nor the vapor- 
liquid transition is a problem, while unmixing occurs for 
[a&lall T <T C = 1.4230±0.005 (unit of T is chosen such 
that e/k B = 1). Working at T = 1.0, the bulk A-rich 
(B-rich) phases are almost pure, x^ ex = (NA/N) coex « 
0.97 (x^ oex = 0.03) 0. 

We choose "antisymmetric" interactions of the walls 
with the fluid particles such that also in the thin film 
geometry phase coexistence between A-rich and B-rich 
phases occurs at chemical potential difference A/i = 
between A- and B-particles, as in the bulk (aH, l39| . 
Specifically, for the A particles we take a wall potential 



u A {z) 



and for the B-particles 
27rp* 



u B (z) 



H(^)+( 



D + 5-z 



where < z < D is the coordinate perpendicular to 
the walls, and the offset 5 = a/2. Both the walls exert 
the same repulsive potential (of strength e r = e/2) on 
both kinds of particles. This potential can be thought 
of as resulting from integration over LJ repulsions from 
atoms residing in a semi-infinite wall. An attractive 
interaction (of variable strength e a ) acts only on the 
A-particles from the wall at z = 0, while the same acts 
only on the B particles from the wall at z = D. 

For direct measurement of 6 from the inclination of 
A — B interface, Monte Carlo (MC) simulations are car- 
ried out in the canonical ensemble where in addition to 
the standard particle displacement moves, exchange be- 
tween randomly chosen particle pairs was also tried. For 
the displacements, a trial shift of a cartesian coordinate 
of a randomly selected particle in the range [— <r/20, 
+a/20] around its old position is chosen and the stan- 
dard Metropolis criterion applied. On the other hand, 
thermodynamic integration (see below) was performed 
by using the raw data obtained from MC simulation in 
the semi-grand-canonical (SGMC) ensemble. In the lat- 
ter, after every 10 displacement steps per particle N/10 
particles are randomly chosen in succession and an at- 
tempted identity switch (A — >> B — >• A) is made (for 
more technical details, see jai|,[a9|). In addition, succes- 
sive umbrella sampling was performed using SGMC to 
obtain inter facial free energies of the co-existing phases. 

As a first step, we estimate jab from the finite- 
size extrapolation of the distribution Pl(xa) (xa — 
Na/N, the mole fraction of A-particles) using a L x 
L x L geometry, with periodic boundary conditions 
in all a>, y- and z- directions. Applying successive 
umbrella sampling [40|, the effective excess free en- 
ergy density /l(xa^T) is found from /l(xa^T) = 
-(l/V)\^^ V Tix A )/P^NVT(x^)l As is well- 
known hjj [aH, |41|, the flat maximum, /hump? in 
Jl(xa,T), near xa = 0.5 in Fig. [TJ is due to the for- 
mation of two planar L x L interfaces in the simulation 
box, separating a slab of A-rich phase from the B-rich 
phase (or vice versa). Thus this barrier is given by 
2jab(L)/L and the extrapolation of such data towards 
L — )> oo yields jab = 0.722 ± 0.002, as demonstrated in 
Fig. Kb). 

Interfaces in nano- confinement: Using the model 
in a L x L x D geometry, applying the wall potentials 
d!]) and (j5j) and starting with initial configurations pre- 
pared such that two A — B interfaces meet the walls 
vertically, we stabilize two-phase configurations of the 
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FIG. 1: (a) Effective free energy density fh{xA,T) of the 
symmetric LJ mixture at T — 1.0 for a box of linear dimen- 
sion L = 10 (in units of a) plotted vs. xa- The estimation of 
the size-dependent interfacial tension jab(L) is indicated, 
(b) Extrapolation of 'Jab(L) as function of 1/L. 

system via MC simulation in canonical ensemble for 
varying strengths of attractive wall potential e a . For 
xa = 0.5 and e a < 0.2 we find slab- like configurations 
of the B-rich phase as seen in Fig. [2j where the inter- 
faces are represented by thick lines. These interfaces 
are perpendicular to the x?/-plane only for e a = when 
the contact angle = tt/2 by symmetry of our model. 
For e a > they are inclined relative to the x?/-plane. 
We extract an effective contact angle # e ff from the av- 
erage slope of these interfaces, measured from the cen- 
tral part of the film, excluding regions of width 2a at 
both walls. For e a = 0.25, wetting has already occurred 
(#eff = 0) and so, instead of a slab with two interfaces, 
we now find a single interface, separating the A-rich do- 
main near the bottom wall that prefers A from 5-rich 
domains near the top wall preferring B. One conse- 
quence of nano-confinement that one can see in all cases 
is a pronounced layering of the total density parallel to 
the confining walls and extending throughout the thin 
film. If one chooses suitable off-critical concentrations, 
e.g., xa — 0.12, one can observe wall-attached droplets 



£a=0.05 




X 



FIG. 2: Snapshot pictures of the configurations of the A- 
particles (dots) in a L x L x D system with L = 32, D — 
8, projected onto the xz-plane, for the case of interfaces 
running along the ^/-direction. B-particles are not shown. 
Four choices of attractive wall-particle interaction strength 
Sa are included, as indicated. The instantaneous interface 
position (averaged over the y coordinate) is highlighted by 
thick lines. 



S a = 




e a = 0.05 




FIG. 3: Three-dimensional snapshots of A— particle droplet 
configurations for the cases xa = 0.12, s a = (D = 12, 
L = 24; upper part) and x a = 0.12, e a = 0.05 (D = 12, 
L = 30; lower part). 
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FIG. 4: Density profiles (pa(z)), (pb(z)), obtained from the 
semi-grand-canonical MC simulation in a LxLxD geometry 
averaging over x and y coordinates, plotted vs. z for a film of 
thickness D — 10. Three top panels (e a — 0, 0.15 and 0.25) 
correspond to the A-rich phases. The panel at the bottom 
with s a = 0.3 shows the profiles for a state where unbinding 
of the interface from the wall at z — D has occurred, with 
the delocalized interface now being at z = D/2. 



shown in Fig. [3l similar to the Ising case 

MM- The 

analysis of such states when obtained via umbrella sam- 
pling method, thus containing the information of /l, 
allows the estimation of the line tension, when the ap- 
propriate contact angle is independently known. 

Thermodynamic Integration and Bulk Con- 
tact Angle: We consider a film of large but finite thick- 
ness in the semi-grand-canonical ensemble with A/i = 
and vary s a being in the regime of incomplete wetting, 
where the system is either in the pure A-rich phase or 
in the pure B-rich phase, prepared via SGMC simu- 
lation. From the density profiles of such systems, as 
shown in Fig. [H where the three upper panels corre- 
spond to A-rich phases (note that in this case only a 
slight enhancement of B particles could be seen at the 
wall preferring the latter), one can extract information 
on j w a — JwB by a thermodynamic integration method. 
To derive this, we first express the total free energy F of 
the film via the partition function as a configurational 



F = -k B T\nZ 

= -k B T\n J dXexp[-/3n b (X)-/3n r w (X) 
+/fe a L 2 (2W3) x 



D 



{pa{z)) (tTs) dz + 


D 



(6) 



Here f3 = l/fc#T and X includes all configurational de- 
grees of freedom of the system. For the sake of brevity, 
neither the bulk energy T-L^X) nor the energy H^X) 
due to the repulsive part of the wall potential has been 
explicitly written down. The density profiles (pa,b(z)) 
are averaged over x and y coordinates. We now take 
a derivative with respect to e ai which singles out the 
derivatives of the wall free energies f s at z = and 
z = D, respectively (f s is normalized per unit area) 



(*=0) , 



D 



B Az=D) 

(irl = (V/3)x 



D 



I d <DTs^) 3{ ^ z)) - (8) 



These equations are readily integrated to yield 
/i z=0) (^) = /i z=0) (0) + (27rp73)x 

s a D 

J K J dz(-^Y <p A (e' a ,z) >, (9) 





/f D) W = r i)) (0) + (2^73)x 

(J 



J de' a J dz( 





D + 5-z 



) <p B (e' a ,z)>. (10) 



Noting that surface free energies of bulk B-rich phases 
are related to those of the A-rich phases simply by sym- 
metry: 

fs ^{^a) \b — rich phase fs ^ (^a) \ A— rich phase? (H) 
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FIG. 5: The full curve is a plot of contact angle vs. e a , 
as obtained from Young's equation with the information 
about j w a — JwB being provided from thermodynamic inte- 
gration, Eq. (f]"2|) . Circles show corresponding observations 
from Fig. El 

we immediately obtain the desired difference 

IwA — JwB = f^ Z ~°\e a ) |A-richphase ~ 

fs H^o) | B— rich phase 
e a D 

= (2np*/3) J de' a j ' dz x 
o o 

(PB(e' a ,z)) A _ iich ( D+ a s _ z ) 3 _. 

(12) 

In Eq. ([T2]) both the profiles {pA{s' a ^z)) and {pB(s' a ^z)) 
are sampled, without loss of generality, in the A-rich 
phase, as shown in Fig. HJ 

Given the knowledge of both the difference j w a — 1wB 
(as function of e a ) and of Jab, we can predict the depen- 
dence of the contact angle on e a , as shown in Fig. [5] 
Following this we obtain the wetting transition (6 = 0) 
at e a = 0.240 ± 0.005. Note also that in our model, for 
e a = we must have = 7r/2 due to the full symmetry 
between A and 5 that is present then. Interestingly, the 
contact angles predicted from Young's equation, which 
refer to quasi-macroscopic droplets, agree almost per- 
fectly with the contact angles, represented by circles in 
Fig. [5j observed in a nanoscopically thin film that is 
only 8 Lennard- Jones diameters thick (Fig. [2}. Thus 
Young's equation still is useful at the nanoscale! 

Estimation of the line tension: It would be wrong 
to conclude from this success of Young's equation that 
the line tension does not play any role in our model. 
In fact for the case e a = it is relatively easy to esti- 
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FIG. 6: Effective interface tension [estimated via 
^/hum P /(2/csT)] plotted vs. 2/7}, for systems of linear 
dimensions L x L x D, with L = 30 and using several 
choices of 7), for the case s a = 0. Arrow shows the esti- 
mate of jab from Fig. [TJb). The straight line is a fit to 
the form jab + 2r/D, where both jab and r have been 
used as adjustable parameters. Fitting gives r = —0.52 and 
jab = 0.725 [note the agreement with the estimate from 
Fig. Kb)]. 



mate it from the flat free energy maximum, /hump (cf. 
Fig.DJa)), for systems in L x L x D thin film geometry 
by varying D. Then one expects that /hump has a value 
2j AB / L + 4t/(LD); hence a plot of Lf hump /(2k B T) vs. 
2/D, as presented in Fig.[6j should yield a straight line, 
where the ^-intercept is the already known interfacial 
tension jab, while the slope yields an estimate of r. In 
this way we find r w —0.52 for s a = 0. Recall that by 
symmetry, = tt/2 for e a — 0, irrespective of the pres- 
ence of the line tension. Unfortunately, when e a > 
and the interfaces are inclined (Fig. [2]), large statistical 
fluctuations prevent us from using; this method. How- 
ever, following Winter et al. j36|, HtJ we have analyzed 
the excess free energy of sphere-cap shaped droplets. As 
will be described elsewhere in detail [42|, the absolute 
magnitude of r, obtained from this analysis, decreases 
rapidly with decreasing contact angle. 

Summary: In this letter, we have described a study 
of the contact angle of nanoscale fluid bridges, for a 
model of a symmetric binary Lennard- Jones mixture, 
for the full range from neutral walls up to the wetting 
transition. Supplying the interfacial tensions entering 
Young's equation, we observe that the quantitative 
agreement of the contact angles thus obtained with the 
direct observation from inclined nanoscopic interfaces 
in a 50 : 50 composition is remarkable. While jab 
and the line tension r are extracted from suitable 
system size dependences, j w a — JwB is estimated via 
a new thermodynamic integration method. Clearly, 
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the symmetric Lennard- Jones mixture is a simple 
model system, but we expect that our methods can be 
generalized to more complex models of real fluids. In 
this way, practically relevant information to guide the 
development of nanofluidic devices will come in reach. 
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